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Abstract 

We calculate numerically the logarithmic contribution to the entanglement entropy of a cylindrical 
region in three spatial dimensions for both, free scalar and Dirac fields. The coefficient is universal and 
proportional to the type c conformal anomaly in agreement with recent analytical predictions. We also 
calculate the mass corrections to the entanglement entropy for scalar and Dirac fields in a disk. These 
apparently unrelated problems make contact through the dimensional reduction procedure valid for 
free fields whereby the entanglement entropy for the cylinder can be calculated as an integral over 
masses of the disk entanglement entropies. Coming from the same numerical evaluation in the lattice, 
each coefficient is cross checked by the other, testing in this way the two results simultaneously. 



1 Introduction 

The entanglement entropy, being a measure of the correlations between two subsystems, depends on 
the geometry of the separating surface, becoming for this reason a quantity with a strong geometric 
character. This has inspired a holographic interpretation within the AdS-CFT correspondence framework 
[Tl[2], which provides a purely geometric way of calculating the entropy of conformal field theories. There 
are several results in the literature where the geometrical properties for different sets, dimensions and 
theories are explored both, from the quantum field theory and holographic approaches. Among these, 
the coefficient s of the logarithmically divergent term in the entanglement entropy ^log = s loge, where 
e is the ultraviolet cut-off, for general conformal field theories in four dimensions has been found to be 
proportional to the type a and c conformal anomaly coefficients [3] 

Here x{9V) is the Euler number of the surface, fc^^ = ^I'jllu'^a'n^p is the second fundamental form, 
with i ^ \,2 are a pair of unit vectors orthogonal to dV , and 7^1, ~ 5^^ — n^^nl, is the induced metric on 
the surface. From ([T]), it can be seen that the sphere and the cylinder are sensitive only to one type of 
anomaly coefficient, a or c, respectively. This is [3] 

^=90' 

for the sphere, while for a cylinder of length L and radius R it is 

For spherical sets, this result was later recalculated, numerically and analytically [H [SJ IHl [Zl IH] , using 
different technics, and recently extended to any theory [51. Instead, the validity of the result for 
the cylinder is more subtle. It has been also studied from the holographic point of view in 10 where s 
is found to be proportional to the c conformal anomaly in four dimensions as in (jS]), but it cannot be 
generalized to higher dimensions where other conformally invariant terms can be added to ([T]). 
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In this paper, we study numerically the logarithmic contributions to the entropy for a cylinder in (3+1) 
dimensions for free massless scalar and Dirac fields. The method we use consists first, in dimensionally 
reducing the problem of the three dimensional cylinder to the one of an infinite set of massive fields living 
in a two dimensional spatial disk. The entanglement entropy is then calculated numerically by the real 
time approach [llj where the reduced density matrix is written in terms of free correlators. This is based, 
on one hand, on the method developed by Sredniki for calculating the entropy for scalars in spherical 
sets [H] where the field is discretized in the radial direction in polar coordinates, and on the other, on 
the work by Peschel, whereby the reduced density matrix can be written in terms of correlators for both, 
bosonic and fermionic, solvable lattice systems |13) . Thus, for the numerical evaluation of the entropy 
of a Dirac field in the disk, we extend the method described in [T^ , originally applicable only to scalar 
fields. 

Then, from the expansion of the entropy in powers of (mi?) valid in the large mR limit, we find the 
coefficient s for the cylinder ([3]). This is directly related to the coefficient of the term {mR)~^. In the same 
expansion, we also identify the coefficient of the linear term mR. This last, has been recently calculated 
analytically in |14| , where the same type of expansion is considered for general smooth geometries and 
massive scalar theories. From the holographic point of view, new contributions to the entanglement 
entropy were also found when considering theories deformed by a relevant operator [15], consistently to 
that reported in [T3|. 

We find, for the cylinder, within a porcentual error ^ 0.15% for scalars and ^ 1% for Dirac fields, 

L 



2A0R ' 
= 240i?' 



and for the disk, within a porcentual error ^ 0.05%, 

AS,n^-'^Rm, (5) 

in agreement with [3] and jl4) . Coming from the same expansion, each coefficient is cross checked by 
the other testing both at the same time. This concordance sums support to the general validity of the 
method used in [3]. We will discuss this issue in more detail, later in the Discussion. 

The paper is organized as follows: in the second Section, we discuss the dimensional reduction pro- 
cedure and show the problem of massless fields in the cylinder can be reduced to the one of massive 
fields in a disk. In the third Section, we discuss the approach we use for the numerical evaluation of the 
entanglement entropy in the disk, for both, scalar and Dirac massive fields. In Section 4 we present our 
results and finally, its interpretation in the Discussion. 



2 The cylinder by dimensional reduction 

For free fields, some universal terms in the entanglement entropy in high dimensions can be obtained via 
dimensional reduction technics from results calculated in lower dimensions [TT]. Let us consider a set in 
three spatial dimensions of the form V = D x X , where X is a line on the first coordinate xi, of length 
L, and I? is a sphere in two dimensions (disk). We are interested in the entropy of V in the limit of large 
L. The direction xi can be compactified by imposing periodic boundary conditions xi = xi+ L, without 
changing the result of the leading extensive term. For a free massless field we Fourier decompose it into 
the corresponding field modes in the compact direction 

<j)nixi,X2,X3,t) = 4>n{x2,X3,t)e'~l^ . (6) 

The problem then reduces to a two dimensional one with an infinite tower of massive fields. For 
example, for a free scalar we obtain the tower of fields 0, 

DiCbn - Ds^n + i^ f^n ■ (7) 
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From the point of view of the non compact X2, directions, these fields have masses given by 

ml=(^-^n\\ (8) 



^ L 

Summing over the contributions of all the decoupled 2 dimensional fields we have 

oo 

S{y)= S{D,m^). (9) 

n— — oo 

In the limit of large L we can convert this sum into integral 

S{V) = - / dmS{D,m). (10) 
^ Jo 

The universal terms in S{V) will then come from the ones of S{D) after integrating over the mass. For 
a Dirac field, the spin multiplicity factor 2[('^+-^)/^l has to be incorporated. In the present case, it is 

S{V) = 2- / dmS{D,m) . (11) 
Ti" Jo 

Expanding the entanglement entropy S{D^m) in powers of niR for large mR 

S{D,m) = cimR + co + c,^)^ + ..., (12) 

^ ' mR 

and inserting (fT2|) in (fTO]) and (fTTj) for scalar and Dirac fields respectively, we obtain for the (3 + 1) 
dimensional theory that the logarithmic coefficient s in S{V) is directly related to C(_i) by 

= -2cf„,A. (14) 

On the other hand, for the dimensionally reduced theory, the contribution in the entropy proportional to 
the mass is given by the term proportional to the coeHicient ci . This contribution was calculated in |14] 
for massive scalar fields in any dimension for a waveguide geometry with specified boundary conditions 
using heat kernel methods. The terms extensive in the area ^d-i hi even spatial dimensions are given by 



12 Je^ t (47rt)- 
From (IT5|) it follows, 

= 7d m''-' ^d-i , (16) 
where = (-l)(''/^Hl2 (27r)('^-2)/2(^_ ^^^nj-i j^^. ^ ^^^^^ Yot Dirac fermions with 2[(''+i)/2] components, 
the same calculation can be done taking into account the extra factor 21^'^+^^/^!^^ [16j, which relates the 
scalar and Dirac 7 coefficients, 

^f^2^i^+im-i^s, (17) 

In our case, d — 2 and 

27r 

ASm = ciR m , ci = 27r72 = ~~[^^ (1^) 

both for Dirac and scalar fields. 



3 The disk: Numerical evaluation for massive scalar and Dirac 
fields 

We evaluate numerically the entanglement entropy for a two dimensional disk. In this section we describe 
the numerical method and models in the lattice for scalar and Dirac massive fields. 



3 



3.1 Massive scalars in a disk 

Consider the quadratic Hamiltonian for a massive scalar in (2 + 1) dimensions 

^=11 dVUMf + {^4>f + m202) . (19) 
Following |12j . we separate variables in polar coordinates and introduce new ones 

k{r) = .[^ fdee"^''mr). (20) 



2tt 

r 

2^ 



^n{r) = I d9e^"\{0,r), (21) 



(22) 



such that 

[Mr),^n'ir')]=iSnn'Sir-r'). (23) 
Then, the Hamiltonian H — J^n ^« takes the form 

1 



For a general quadratic discrete Hamiltonian H = ^^T^f + \ X^ij 4'iKij(f>j, (dropping the n index 
temporarily), the vacuum (ground state) two point correlators X^j and 

(M,) - X,,, (25) 

{7^^7T,) = Py, (26) 

are given by [H] 



X,, = {^,<j>,) ^ -{K-^),, , (27) 



Py = (7r,^,) = -(i^5),,, (28) 

in terms of the K matrix. Here, I < i, j < N , where N is the size of the lattice and acts as an infrared 
regulator. 

Then, the entropy of the disk is given by 

oo 

SD =So + Y^ 2Sn , (29) 



n=l 



with [TT] 



5„ ^ tr [{.JX^PJ^ + -) \ogi^ X^PD + _) _ i^X^^Pj^ - -)\og{^ XDPD _ _) j . (30) 

The superscript D in (15(11) means the indices of the matrices are restricted to the region D: I < i, j < r, 
with r the radius of the disk. 

In this case, after discretization we find the matrix Kn corresponding to is 

= l+n^ + m\ (31) 

9 

K = 2 + -^+m^ (32) 
i^M+i = - ' + ^/^ =K:+^'\ (33) 

V7(7TT) 

Summarizing, the numerical evaluation of the entropy for massive scalar fields in a disk starts with the 
calculation of the {N x N) matrix Kn for a given mass m, angular momentum n, and infrared regulator 
A'' which gives the size of the unidimensional lattice. From (|27|) and psp we calculate the two point 
correlators X and P. Then, we reduce them to the disk and calculate the contribution Sn (|30|) . Finally, 
the entropy is given by the sum 
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3.2 Massive fermions in a disk. 

The Hamiltonian for a massive Dirac field in (2 + 1) dimensions can be written as 

H = j dVi>\x,y)H^^^{x,y) , (34) 
with the one particle Hamiltonian 



H^^^ =i^^-{aA + avdy)+^m. (35) 



Choosing 



Q-x = CTi , (36) 

(Xy = (72, (37) 

/3 = C73, (38) 

with Gi the Pauh matrices, the relations {ai,aj} = {aj,/3} = and a? = = 1 are satisfied. In polar 
coordinates r, 6 the Hamiltonian takes the form 



i?^^^ = T \-hede + hrdA + mhm , 
I \r j 



(39) 



where he, hr and hm are 



id 



hr = ( ^te \ ) , (41) 
hm = { i ^,]. (42) 



-1 ^ 

Separating variables, we propose the following two component Dirac spinor 



where (/)^ — -^|=e'^^"+5) and (/)^ = -^e'^^""^) are the eigenvector components of the angular momentum 



operator J = ide — has with half integer eigenvalue n 



Note that J commutes with H. Then, we can express the Hamiltonian as a sum over n such that 

H = Y,Hn = Y. f drr ai{r)H^an{r) (45) 

where 

H^^^ = — (iai/2 + nc72) - iaidr + mcia , (46) 



r 

and 



""(0 = ( !"S ) • (47) 



^'r^(^') 

In general, in the discrete free case, the reduced density matrix can be written in terms of a Hermitian 
operator H 

p = Ke-'^ = Ke^v ^l^i^^i . (48) 
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This can be expressed in terms of the correlators Cij 

(^.vj) = aj , (49) 

via the identification [11] 

•H = -log(C~i - 1). (50) 

Here, the tpi a-re fermion operators canonically normalized. Then, the entropy can be written in terms of 
C as = -tr((l - C) log(l - C) + ClogC). 

For a general quadratic case, with discrete Hamiltonian H — Mijipjipj, the correlator is directly 
related to My by 

C = e(-Af ) . (51) 

In (|5ip . the indices run over the complete space. Then, we introduce an infrared regulator N which is 
the size of the lattice. 

Therefore, the entropy S for the disk D, is given by a sum over the angular momentum n as 

S{D) = J2 -tr((l - C,?) log(l - ) + logC,f ) . (52) 

n 

The correlators for fixed angular momentum n are restricted to the disk region D. 

In order to identify the M matrix in the case of the disk, we first introduce the operators a(r) = 
r^^/^''a(r). These are normalized such that they satisfy the canonical anticonmutation relations required 
for the application of ([5^ . The discrete Hamiltonian in these variables for a N lattice, takes the form 
(dropping the angular momentum index temporarily) 

N N / _ ^ 

^ = E^'^^-= E("*'^*)M,J , (53) 

i,j=l i,j=l \ 3 / 

where the indices i,j are the discrete variables corresponding to the continuum radial coordinate r. Mij, 
for fixed i, j, is a (2 x 2) matrix such that H,^j = Mlju*Uj + Mlfu*Vj + Mfji*Uj + Mffi*ij with 
i,j = l,...,iV. 

Finally, we define a 27V x 27V matrix Mlk+a-2,2i+fi-2 ^ ^ f^j. fc^ ; ^ iV and a, /3 = 1, 2. The 
non zero entries of M (for each angular momentum n) are 

Mt = (-1)'=+V, (54) 
Ml'^ = ^{n+\) , M2.i = -z(„+i), (55) 

^2fe~l,2fe-2^^ , M2Mfe+l^_i^ (5g) 

The Mn matrix satisfies Af,J = M„ and it has real eigenvalues symmetric with respect to the origin. 
There is also a symmetry related to the mapping n — —n which can be seen also in the continuum limit. 

Summarizing, the numerical evaluation of the entropy for massive Dirac fields in a disk starts with 
the calculation of the (2iV x 2N) matrix M„ for a given mass m, angular momentum n, and iV, the size 
of the unidimensional lattice. From (fSTj) we calculate the two point correlator C (|49p . Then, we reduce 
it to the disk and calculate the contribution Sn in ([5^ . Finally, the entropy is given by the sum 



4 Results 

We calculate the entanglement entropy for a disk in a lattice of 200 points. We consider regions of radii 
r going from 50 to 80 in lattice units for the scalar and from 30 to 50 for the fermionic field. The sum 
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0.00 0.05 0.10 0.15 0.20 0.25 0.30 

m 



Figure 1: Scalar field: The points correspond to the coefficient of the linear term in r in the disk 
entanglement entropy for different masses. The linear coefficient in m in the fit, drawn with a solid line, 
is -0.52359 f|. This is the value of Ci in 




Figure 2: Scalar field: c^i{m) — — c_i(to). The points correspond to the coefficient of the term ^ in the 
disk entanglement entropy for different masses. The coefficient of the term proportional to ^ in the fit 
drawn with a solid line is 0.01342 ^ This is the value of — c_i in ([T^ . 



2.7 - 



2.4 - 




0.0 0.1 0.2 0.3 0.4 

m 



Figure 3: Dirac field: The points correspond to the coefficient of the linear term in r in the disk entan- 
glement entropy for different masses. The linear coefficient in m in the fit, drawn with a solid line, is 

-1.04658 2 X 2|. This is ci in (HH). The extra factor of two is due to the fermion doubling in the 

lattice. 
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0.1 0.2 0.3 0.4 0.5 

m 



Figure 4: Dirac field: c_i(m) = — c_i(to). The points correspond to the coefficient of the term i in the 
disk entanglement entropy for different masses. The coefficient of the term proportional to ^ in the fit, 
drawn with a solid line, is 0.07754 ^ This is the value of — c_i in (fT2|) . The factors 2 and 1/2 due 
to the fermion doubling and the spin multiplicity in the dimensional reduction formula (jl4p respectively, 
compensate each other. 

over the angular momentum n is done exactly up to rimax — 3000 and the corrections coming from 
contributions n > rimax are added by fitting the exact entropy for some large values of n. The fit we use 
is of the form 

1 , logn 1 logn 1 logn 

Sn = 0,2^ + 02 5- +04^+04 J- + a^—^ + bg tt- (59) 

rt^ n° n° 

since for large n {n N) both C and \J XP defined in ([51]), ([ST)) and ([28]) have expansions in even 
inverse powers of n [12]. 

This is done for different masses in the range 1/20 <m< 1/2 such that mT^ > e (we fix e = 1) and 
mr > 1. For each mass, we fit the entropy in terms of the size r of the disk, 

S* = ci(m)r + co(m) + c_i(m) — h ... . (60) 

r 

Once the coefficients ci(m) and c_i(m) are identified, we expand in powers of m, 

ci(to) = cim + cj + Cj^^— , (61) 

c_i(to) = cLim + c° 1 + c_i^ . (62) 

We obtain ci and c_i of as the m and ^ coefficients in (IFTI) and (15^ respectively. For example, 
in the scalar case, for m = 1/10 we calculate the entanglement entropy for disks with radii from 60 to 
80 (in lattice units). From the best fit, we extract the coefficients of the terms proportional to r and 
1/r. In this example, this gives 0.20744 - 0.13303/r + 0.41493r. Then, ci(m = 1/10) = 0.41493 and 
c_i(m = 1/10) = —0.13303. The values of ci(m) and c-i(m) are shown in Figures (1),(2) for scalars and 
Figures (3), (4) for the Dirac field. In order to compare the lattice results with the continuum expectations 
for fermions, a factor 1/2 has to be incorporated in the entropy formula due to the fermion doubling in 
the unidimensional radial lattice. We obtain 

7r 

c^ -^^ — —0.01342 ^ —^^ , within a porcentual error of 0.15 , (63) 



0.07754 ^ , within a porcentual error of 1.2 , (64) 



(-1) ■ 240 



and 



27r 

-0.52359 ^ ' within a porcentual error of 0.02 , (65) 
27r 

-0.52329 ^ 1 within a porcentual error of 0.05 , (66) 
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in agreement to the analytical results. 



5 Discussion 

Our numerical results agree with the ones predicted analytically within porcentual errors 1.2% or less. 
Since both results, the one for the cylinder and the one for the disk, come from the identification of 
different coefficients in the same expansion (fT^ . they cross check each other giving a solid basement to 
our numerical test. 

On the other hand, by the analytical side within the QFT approach, the general grounds the result 
for the cylinder is based on, rely first on the application of the replication method, then on demanding 
conformal invariance of the logarithmic contribution in the entropy derived from the effective action of 
the replicated manifold, and finally on calibrating a free parameter by using the holographic correspon- 
dence [1]. The replica trick is a standard approach to calculate the entanglement entropy [17]. In this 
construction, the trace of the integers powers of the reduced density matrix corresponds to the partition 
function of n copies of the system connected consistently through the cuts along the boundary. The 
replicated manifold is non trivial due to the conical singularities of deficit angle a = 27r(n — 1) placed on 
the boundary. From this, the entropy associated to the set V can be calculated as 

S = -lim ^tr[p'^>] ^ - hm (— - 1) log Z„ . (67) 

n^l On n^l On 

In the above formula we assume that not only the partition function Z„ is calculable, but it is also possible 
to analytically continue it to non integer n. The method presented in [3] differs subtly from the above: 
instead of continuing analytically in n the partition function Z„, they first analytically continue the cover 
manifold assuming the corresponding geometry can be defined for an arbitrary and small angular deficit 
in the limit n — )■ 1 and then calculate the corresponding Z„ [3l [18]. Taking into account this is not 
possible in general, then, the identification of the logarithmic coefficient with linear combinations of the 
stress tensor anomalies is not fully justified [lOl [191 [20l [21] . In the case of spherical sets, the objections 
to the applicability of the replication method disappear since there is an extra rotational symmetry in 
the transverse space about the entangling surface. In fact, the result for the sphere in four dimensions 
can be extended to any dimensions and any field theory [9^ . These results complete a solid proof for the 
term proportional to the type a anomaly in ([T]). 

On the other hand, for the cylinder, the situation is more obscure. Very recently, it was found in |10] 
for holographic models, that in four dimensions the logarithmic coefficient is proportional to the anomaly 
as expected, but that the same is not true for higher dimensions, where new conformally invariant terms 
could be added to ([1]). In this scenario, we conclude our results give support to the validity of the 
Solodukhin's calculation for the cylinder in four dimensions. 
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